An Approach for Delaunay Tetrahedralization of Bodies with Curved Boundaries

Problem of tetrahedral meshing of three-dimensional domains whose boundaries are curved surfaces is wide open. Traditional approach consists in an approximation of curved boundaries by piecewise linear boundaries before mesh generation. As the result mesh quality may deteriorate. This paper presents a technique for Delaunay-based tetrahedralization in which a set of constrained facets is formed dynamically during face recovery and mechanisms for mutual retriangulation of the curved faces and the tetrahedralization are suggested. The proposed algorithm is constructed in such a way that a facet that was once added in the set of constrained facets is never split into small triangles. It allows retaining the high quality of surface mesh in the tetrahedralization, because during boundary recovery the surface mesh on the curved faces and the tetrahedralization are refined conjointly.